Magnetic edge states of impenetrable stripe 
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The electron motion in a strong perpendicular magnetic field close to the impenetrable stripe 
is considered by making use of the singular integral equations technique. The energy spectrum is 
calculated and compared with the energy spectrum of the round antidot. It is shown that in the case 
of the long stripe the eigenfunctions can be obtained as a superposition of magnetic edge modes, 
while fractal energy levels obtained in a high energy region can be explained from the quasi-classical 
point of view. 

PACS numbers: 73.20.Dx, 85.30.Vw, 03.65.-w 



-a 
c 

o 
o 



o 

(N 



(N 
O 

-i— > 
c3 



i 

C 

O 

o 



X 



I. INTRODUCTION 

Progress in a nanometer technology and the ability to 
tailor potentials has triggered a broad activity in low- 
dimensional semiconductor nanostructures. Among the 
2D (two-dimensional) structures the quantum dots with 
the electrons confined in a small region have been a sub- 
ject of intense theoretical and experimental research dur- 
ing last years 111 The complete confinement and the dis- 
crete energy spectrum converted these objects into a use- 
ful instrument for the electron interaction and correla- 
tion studiesfl In the strong perpendicular magnetic field 
the quantum antidot, the region with a repulsive poten- 
tial, can bound the electrons as well. The magnetotrans-j 
port experiments on the arrays of the quantum antidotsS 
showed the close relation of the pronounced structure in 
the magnetoresistance and the periodic classical orbits, 
or the corresponding spectrum of the antidots. It was 
tested on the arrays of various shape antidotsrl The im- 
portance of the antidot-bounded electron states was con- 
firmed in the studies of magnetotransport through clus- 
ters of the antidotsa and the individual antidots E 

The spectrum of the quantum antidot in the magnetic 
field is also interesting from the dynamic chaos point 
of view (see review articldl and references there). The 
quantum antidots together with quantum billiards are 
the most simple and convenient structures for revealing 
the links between the auto-correlations in quantum spec- 
trum and the periodic orbits of the classical problem. 
The most convenient technique for solving the antidot 
eigenvalue problems with not separating variables is the 
singular integral equations. Usually the sharp antidot 
edges increase the singularity of these equations making 
them rather complicated and even not useful. In this 
paper using the simplest antidot, the finite impenetrable 
line, we demonstrate how the integral equation technique 
can be used in the case of the antidot with the sharp 
edges. We have failed to find that this simple but reveal- 
ing antidot spectrum has been considered ever before. 
Comparing this spectrum with the round antidot one we 
demonstrate the main features of the problem with non- 
separable variables. The limit cases of a long antidot 
and high electron energy show the peculiarities of the 
quantized magnetic edge modes and the quasi-classical 
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FIG. 1: Layout. 



quantization. 

The paper is organized as follows. In Sec. II the prob- 
lem is formulated. In Sec. Ill the simplest ultra short 
stripe case is considered, and in Sec. IV the numerical 
results for the energy spectrum are presented. In order 
to explain the physical meaning of the spectrum pecu- 
liarities two limit cases — the long stripe in Sec. V and 
large electron energy in Sec. VI — are considered. In 
the last Sec. VII the conclusions are given, and in the 
Appendix the details related to the discretization of the 
singular integral equation are collected. 



II. MODEL 

We consider the electron moving in iry-plane which is 
shown in Fig. 0, The antidot, an infinitely thin impen- 
etrable line \x\ ^ a/2, y — 0, is indicated by the grey 
stripe. We solve the Schrodinger equation 

{H - E}^{r) = (1) 

with the following dimensionless Hamiltonian: 



if = -i{V + ? :A(r)} ; 



(2) 



where the perpendicular magnetic field is described by 
the vector potential in the symmetric gauge A(r) = 
{— y,x}/2. We use the following notation for 2D vectors 
r = {x, y}. The energy is measured in hu c (uj c — eB/mc) 
units, and the coordinates — in the magnetic length 
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Ib = y/ch/eB units. The antidot — the impenetrable 
stripe — is taken into account by hard wall boundary 
condition 



$(r) = 0. 

\x\^a/2,y=±0 



(3) 



Besides, the wave function satisfies zero boundary condi- 
tions at the infinity ^(r^—nx, = 0. 

Making use of the Green theorem the above two- 
dimensional problem can be transformed into one- 
dimensional integral equation. Indeed, introducing the 
whole plane Green function as a solution of the equation 



{H - E}G{v\v') = -8{r-r') 



(4) 



and taking the boundary conditions G(oo|r') = into 
account, one can present the wave function as the integral 
(see the details in Ref. [?]) 



1 /-a/2 

tf(r) = - / dx'G(r\x',0)F(x') 

2 J -a/2 



over the perimeter of the antidot. Here 

F{x) = V v (x,+0) - %{x,-0) 



(5) 



(6) 



is the difference of the wave function derivatives on the 
opposite sides of the stripe. We shall refer to it as the 
perimeter function. 

The wave function defined via Eq. (||) satisfies already 
equation (Q) and the boundary condition at the infinity. 
Satisfying boundary condition on the stripe (§) we get 
the following integral equation: 



a/2 



a/2 



dx'K(x,x')F(x') = 



with the kernel 



K(x,x') = 27rG(x,0|x',0). 



(7) 



(8) 



The nonessential factor 2ir is included for the sake of 
convenience. 

This integral equation is our main instrument. Note 
that taking the derivative of Eq. (||) over y and equating it 
to the perimeter function on the stripe, one more integral 
equation can be obtained. In our case it isn't necessary, 
because our stripe-antidot has no inner region, and con- 
sequently, there are no spurious eigenstates, which have 
to be properly eliminated in the case of other antidots.H 

In order to fix the kernel we have to solve Green func- 
tion equation (^) . The solution of it is knownQ 

G(r|r') = --Lexp{i[r xr'] 2 /2} S ( S ), (9a) 

Z7T 

g(s) = -r(-e)exp(-»/2)tf(-e|l|a), (9b) 

J\2 



|r- 
E 



vf/2, 
1/2. 



(9c) 
(9d) 



Here the symbol T(z) stands for T-function, and C/(a|6|z) 
is the Kummer function of the second kinds- the solution 
of the confluent hypergeometric equation.u 

Inserting the above expression into Eq. (||) we obtain 
the following kernel 



K(x,x')=-g(s), 



„'\2 



(10) 



III. SHORT STRIPE 

When solving integral equation (Q) numerically the 
main problem is the kernel singularity at x = x' . One can 
expect that it will lead to the singularity of the perime- 
ter function at the ends of the stripe x = ±a/2. In order 
to reveal the above singularity we considered the limit 
case of the ultra short stripe (a — * 0) when the Kummer 
function can be replaced by its following expansion: 



{Ins + </>(-£) -2V>(1)}. (11) 



Here the symbol if){z) stands for the logarithmic T- 
function derivative, or the so called ^-function. Note 
apart the singular logarithmic term we included terms 
which are large close to the first Landau level (e -> 0), 
and the constant term ip(l) = —7 (7 ~ 0.5772), which 
we need to get the proper behavior of the lowest antidot 
energy branch. 

Now inserting the above expansion into Eqs. (|Io|), then 
into (0), and scaling the variables x — ► ax, we arrive at 
the following approximate integral equation for the ultra 
short stripe case: 

1/2 ,1/2 

dx'\n\x-x'\F{x') = \ I F(x')dx', (12) 
-1/2 J-1/2 



where 



A = --{ln( a 2 /2)+VX-£) + 2 7 } 



(13) 



Thus, the short stripe energy spectrum problem is re- 
duced to the calculation of the eigenvalues of integral 
equation ([l2] ) . Then the energy can be obtained by solv- 
ing algebraic equation (|l3|). For instance, replacing the 
i/'-function by its simplest expansion in the vicinity of the 
first Landau level VK~ e ) ~ V e — 7 we S e t 



E- 1/2 = s 



7- 2A- ln(a 2 /2)' 



(14) 



Due to the average of the eigenfunction over the stripe 
on the right-hand side of equation (12) it has the single 
not equal to zero eigenvalue A. It can be checked by the 
straightforward integration that the function 



F(x) = l/Vl/4 - x 



(15) 



satisfies integral equation (12) with A = —2 In 2. The 
corresponding energy spectrum branch is shown in Fig. 



FIG. 2: Energy spectrum of ultra short stripe. Solid curve - 
solution of Eq. (]l 3|), dotted curve - simplified asymptotics ac- 
cording to Eq. (p"3|), and dashed curve - numerical integration 
of integral equation (Ym. 
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| as a function of the stripe length a: the solid curve 
is obtained by the numerical solution of Eq. ( |l^ ) , while 
the dotted curve indicates the simplified version of the 
asymptotics according to Eq. Both of them co- 

incide in the limit case a — ► 0. We see that the anti- 
dot (short impenetrable stripe) expels a single level from 
the first degenerate Landau state (indicated by a thick 
dashed horizontal line). The longer the stripe is, the 
higher is the level. Note a rather fast energy grow at 
small a values. It is a characteristic feature of the energy 
level corresponding to the non-perturbed level with zero 
orbital momentum, which has a non zero electron density 
at the origin r = 0. All other non-perturbed levels have 
zero electron density there, and thus, they are weakly in- 
fluenced by the stripe, and consequently, not expelled in 
this simplest ultra short stripe approach. 



IV. NUMERICAL RESULTS 

The most important result for the ultra short stripe 
case presented in the previous Section is perimeter func- 
tion ( |T5| ) which singularity at the stripe ends is caused 
by the interplay of the logarithmic singularity of the ker- 
nel and the sharp antidot edges. Thus, it is inherent to 
the perimeter function of general integral equation (Q) as 
well. That is why in order to achieve the proper accu- 
racy in numerical solution of the above equation one has 
to take both singularities (kernel and perimeter function) 
into account explicitly. For this purpose we have replaced 
the perimeter function as follows: 



F(x) = 



fix) 



(16) 



FIG. 3: Energy spectrum of the stripe. The non-perturbed 
Landau levels are indicated by thick dashed horizontal lines. 
The solutions of Eqs. (^) and ( f38| ) are indicated by dotted 
curves. 



proper weights of the discretization scheme. See the de- 
tails of the calculation in Appendix [A|. Instead of solving 
the obtained matrix equation 



ACf = 



the corresponding eigenvalue problem 



(17) 



(18) 



and discretized the obtained integral equation for the 
function f(x) including the singular factors into the 



for various electron energies e was considered. The elec- 
tron energy was defined by zeroing the obtained eigen- 
values A„ = A„(e) = 0. 

The obtained stripe energy spectrum is shown in Fig. || 
where six highest levels expelled from each Landau level 
(dashed horizontal lines) are indicated. On the axes the 
original dimensions are shown. Thus, not only the energy 
dependence on the stripe length but its dependence on 
the magnetic field strength (a/Is ~ VB) can be traced 
as well. 

For a comparison the spectrum of round impenetrable 
antidot with a diameter a is shown in Fig. |. It was 
obtained by means of zeroing the radial antidot wave 
function (it coincides with Green function (^|) with r' = 
assumed) at the antidot border r — a/2. 

The characteristic feature of both spectra is fast ex- 
pelled first antidot level for each Landau state in the case 
of small a values. The detailed behavior of this level ex- 
pelled from the first Landau level for the stripe is shown 
in Fig. H by a dashed curve. Rather good coincidence of 
it with the short stripe energy (solid curve) confirms a 
good accuracy of the developed numerical scheme. 
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FIG. 4: Energy spectrum of the round antidot. 



We see that these two spectra of the antidot-stripe and 
the round antidot differ essentially. The round antidot 
spectrum is a typical one for the system with separable 
variables. In this case the variables can be separated 
due to the cylindric symmetry of the problem, and actu- 
ally we have independent radial problems for every an- 
gular momentum value, which energy spectrum branches 
freely intersect each other. The main point is that when 
the antidot level with some orbital momentum I reaches 
the next Landau level, the antidot level with the same 
momentum I is already expelled from it. Consequently, 
any antidot level freely crosses any Landau level. And 
we see the energy spectrum branches going up when a 
increases with numerous crossings. 

This is not the case for a stripe spectrum. Due to 
the lack of symmetry the orbital momentum is not good 
quantum number any more, and instead of crossings we 
have anti-crossings. Moreover, expelled antidot levels can 
not cross Landau levels any more. So, when the parame- 
ter a/Is grows the expelled antidot levels saturate bellow 
the next Landau level. 

Nevertheless there are still some crossings. See, for in- 
stance, the behavior of levels expelled from the second 
Landau level in Fig. ||. The matter is that the stripe in 
the perpendicular magnetic field still has the inversion 
symmetry (r — > — r). Due to it all the perimeter func- 
tions (and the wave functions as well) can be divided into 
the symmetric and the anti-symmetric ones which actu- 
ally satisfy the different integral equations. Thus, in the 
energy spectrum of the stripe crossings between the sym- 
metric and anti-symmetric spectrum branches are possi- 
ble. These crossings and the occurring waviness of the 
spectrum branches are the most prominent feature of the 
considered stripe spectrum. Note that when the electron 
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FIG. 5: Spectrum of the magnetic edge modes. 

energy grows the above waviness is transformed into nu- 
merous pronounced plateaux on the spectrum branches. 
Now we are going to explain the physical meaning of 
these plateaux and waviness considering two limit cases 
of the long stripe and large electron energies. 



V. LONG STRIPE 

Let us start with the long stripe approximation. When 
the stripe is infinite there are two types of magnetic edge 
modes propagating along the stripe on both its sides. 
In this case we have the problem invariant under the 
translation along the stripe, and thus, the correspond- 
ing eigenvalues can be labelled by the electron momen- 
tum component k along ^-direction. The eigenfunctions 
can be expressed in terms of parabolic cylipder function 
as #(fc|r) = cxp(±ikx)D e {±V2(y - k)}H The spec- 
trum is obtained zeroing the wave function on the stripe 
D e (^fy/2k) = 0, and it is shown in Fig. |[ We see that 
it consists of two independent systems of branches de- 
scribing the electron motion to the left above the stripe, 
and to the right — below it. In the asymptotic region 
k —>■ ±oo the branches tend to the Landau levels shown 
by thick dashed horizontal lines. The intersection points 
of the above branches coincide with the Landau levels as 
well. 

In the case of the finite stripe the electron motions 
above and under the stripe are no more independent, 
because moving above the stripe the electron reaches its 
end, bends around the corner, continues its motion under 
the stripe, and so on. Bending of the corner is a rather 
complicated diffraction problem, but in the asymptotic 
long stripe case we can replace it by some scattering ma- 
trix acting on the longitudinal motion exponents. 

The description of electron motion depends on the 
number of edge modes participating in it. For instance, 
if the electron energy is lower than the second Landau 
level with energy 3/2 (see, the lower thin dotted hori- 
zontal line in Fig. ^ labelled by Eq) there are only two 
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edge modes indicated by solid circles. One of them with 
momentum k moves above the stripe to the left, while 
the other one with momentum — fc moves to the right 
under it. At the stripe ends these edge modes are scat- 
tered one into another. As there is a single scattering 
channel only, the scattering probability is equal to unity. 
Consequently, due to the scattering event the electron 
wave function is multiplied by some scattering ampli- 
tude S = exp{ix(A£')}, while the propagation along 
the stripe can be taken into account by the propagator 
exp(ifca). Thus, taking into account the periodic motion 
of the electron (after bending both stripe ends the expo- 
nential part of the electron wave function must coincide 
with itself), we can write down the following simple rule 
for quantization of the edge modes in the asymptotic long 
stripe case: exp{2i(fca + x)} = 1, or 



ka + x(A£) = 7m, n = l,2, 



(19) 



Now using the relation AE — vk (v — 5/4-^/2) which 
follows from the properties of the parabolic cylinder func- 
tions close to the intersection point at E = 3/2, k = 0, 
and the expansion 



X (AE) = xo + Xi AS, 
we solve Eq. (|l9| ) and get 

A(n - A) 



AE n = 



a + d 



A 



3.55. 



(20) 



(21) 



Two other parameters: the effective elongation of the 
stripe 5 — v\i and the quantum number defect A = xo/ 71 " 
depend on the scattering amplitude phase and unfortu- 
nately, can not be found analytically. 

Fitting the numerically obtained energy branches ex- 
pelled from the first Landau level in the interval 50 < a < 
100 by E = A n /(a + S n ) we have obtained the following 
parameters: 



n 


1 


2 


3 


4 


5 


6 


Sn 


4.29 


5.53 


5.69 


7.25 


7.99 


9.49 


A n 


3.48 


7.04 


10.40 


14.01 


17.39 


20.99 


An/A 


0.99 


1.98 


2.92 


3.93 


4.89 


5.88 



Note the numbers in the last row coincide rather well 
with the integers n what convince us that the picture of 
quantized edge modes is quite adequate. 

To explain the behavior of energy branches expelled 
from the upper Landau levels is more complicated be- 
cause there are more edge modes present. For instance, 
close to the third Landau level (see the upper thin dot- 
ted horizontal line in Fig. || labelled by E\) there are 
four edge modes. Two of them with momenta —k and 
k' propagate above the stripe to the left, while the other 
two with momenta k and —k' propagate under it to the 
right. Consequently, in this case the propagation of elec- 
tron on both sides of the stripe has to be described by 



the following propagator: 



^ — ikoi 



acting on the state vector 



* = 



(22) 



(23) 



composed of edge mode superposition coefficients. Bend- 
ing of the edges is characterized by some 2x2 scattering 
matrix S. Now the quantization is performed by the fol- 
lowing self-consistency condition: 



SV(a)SV(a)y = 



(24) 



and the energy spectrum of the stripe can be defined 
zeroing the determinant of the above equation 



det \SV{a)SV{a) - 1\ = 0. 



(25) 



The absolute values of the scattering matrix S elements 
are given in Ref. [l(| We see that close to the third Lan- 
dau level the absolute value of the off-diagonal elements 
are nearly unity (|5oi| = \Sw\ « 1), while the diagonal 
elements are small \Sqq\ = \Sn\ « Afc/2fco where the 
symbol fco stands for the edge mode intersection point 
and Ak = k' — k is the difference of edge mode momenta, 
propagating on both sides of the stripe. Thus, adding 
the phases we construct the following scattering matrix 



S 



& / e^Afc/2fc 



iifi 



-e- lv Ak/2k Q e l 



(26) 



which satisfies the unitarity condition SS + — X with the 
accuracy of Afc terms. 

Now expanding the momenta close to the intersection 
points as 



k 

k! -k 



k - —AE, k' = fc + —AE, 

Vi v 2 



-AE, 



-(- + - 

2 I V! v 2 



(27a) 
(27b) 



we transform Eq. (^H|) into the following approximate 
equation 

ex.p(iaAE/v) — exp(— 2i$) = sin(fc a — V>)- (28) 

k v' 

If one neglects the small term on the right-hand side of 
this equation one gets the result similar to Eq. ( ^l| ) 



(0) _ A'(n-X) 
a + o' 



(29) 



where $(AS) = $ + $ X AE, A' = nv, A' = $ /vr and 
5' = <f>i«. The right-hand side term taken into account 



FIG. 6: Electron density for long stripe (a/Is = 48): (a) - for 
the highest antidot level expelled from the first Landau level, 
and (b) - for the same level expelled from the second Landau 
level. 
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as a perturbation leads to the following oscillating cor- 
rection: 



A£W = ± 



AE, 



(0) 



■ sin(fcoa — tp). 



(30) 



It is remarkable that the period of oscillations depends 
on fco = l/\/2 only. It leads to Aa — 2-rt/kQ w 8.89, what 
coincides rather well with the period value 8.8 obtained 
from the numerical calculation result for the two upper 
levels expelled from the second Landau level. 

This simple asymptotic picture of interfering magnetic 
edge modes is confirmed by Fig. ^ where the contour 
plots of the electron densities corresponding to the above 
considered antidot states are shown. We see that the 
wave function for the antidot level expelled from the first 
Landau level looks like a cigar (and it does not matter 
how long is the stripe) what indicates that it is composed 
of a single magnetic edge mode. In the case of the antidot 
level expelled from the second Landau level there are the 
lumps (the longer the stripe is the more lumps there are) 
which are caused by the interference of two pairs of edge 
modes propagating on both sides of the stripe. 



VI. QUASI-CLASSICAL LIMIT 

As it has already been mentioned in Sec. |^ (see Fig. ||) 
the oscillations surveyed close to the lowest Landau levels 
change into well pronounced plateaux when the number 
of Landau level is incremented. The detailed view for 
the antidot levels expelled from the fifth Landau level is 
shown in Fig. 0. It is remarkable that the energy of these 
plateaux is very close to simple fractions of cyclotron en- 
ergy. Thus, the plateau indicated by a dashed horizontal 
line is right in the middle between two adjacent Landau 
levels with energy £ = 9/2 + 1/2 = 5, while the energy of 
two other plateaux indicated by dotted horizontal lines 



FIG. 7: The stripe energy levels expelled from the fifth Lan- 
dau level. 



exceed the fifth Landau level energy (9/2) by one and 
two thirds. 

Unfortunately, in the case of high Landau levels the 
number of the interfering edge modes is large, and this 
fact makes it difficult to apply the long stripe approxi- 
mation considered in the previous section. Nevertheless 
the simplified description is still possible due to the large 
electron energy. 

It is known that when the electron energy is large the 
quasi-classical approach based on Bohr quantization rule 
can be used. In the case of free 2D electron in the ho- 
mogeneous perpendicular magnetic field it reduces to the 
estimation of the following integral: 



£dr) = 27m + 7r, n = 0, 1, 2, • 

composed of fast coordinates^ 

£=Px-y/2, V = Py + x / 2 



(31) 



(32) 



over the electron trajectory. Insert ing the solution £ = 
\ // 2Ecost 1 r] = \J2E sini into Eq. ( |3l] ) one immediately 
gets a well known expression for Landau level energy 
E n =n+l/2. Note the considered electron motion is 
two-dimensional, and consequently, two more coordinates 
— the slow motion coordinates 



X = x/2- Py , Y = y/2 + Px 



(33) 



— have to be taken into account. In free electron case 
it is trivial because the Hamiltonian does not depend 
on them. The single important thing is the commutator 
[X, Y] = i, which shows that the total slow coordinate 
phase volume divided by 2ir gives the degeneracy of the 
corresponding energy level. 



FIG. 8: Two circle trajectory. 



FIG. 9: Electron density for the level indicated by short dash 
in Fig. m corresponding to the classical two-circle trajectory. 



The influence of the stripe on the classical electron tra- 
jectory can be taken into account via scattering events. 
The simplest trajectory with two scattering events is 
shown in Fig. || by a solid curve. We see that in this 
case the trajectory is composed of two Larmor circles, 
and consequently, it is twice longer than the trajectory 
of the free electron rotating in the magnetic field. Thus, 
the integral on the left-hand side of Eq. ( |3ll ) becomes 
twice larger. This fact leads to the twice smaller separa- 
tion of energy levels (AE = 1/2) as compared with the 
separation of Landau levels. 

In order to decide whether the fractal Landau levels ob- 
tained in this quasi-classical way can take place or not, 
one has to inspect the slow motion coordinates. We shall 
take them into account in the most simple way. Note that 
there are more equivalent trajectories with two scatter- 
ing events and the same energy (the same radius of the 
Larmor circle) . Let us mark them by the position of the 
left circle center (X, Y). So, changing the center by some 
(AX, AY) we obtain another equivalent trajectory, as it 
is shown in Fig. || by a dashed curve. Thus, the integral 
over all possible coordinates X and Y 



V 2 (R,a) = 




dXdY 



= J dYn 2 (R,2\/R 2 -Y 2 ), (34a) 

n 2 {R, b) = (26 - a)0(2b - a)0(a - b) 

+ aO(b-a) (34b) 

gives the total phase volume for the trajectories with 
given radius R. Here b is the distance between the 
centers of both circles. We assume that the quantity 
N 2 (E,a) — V 2 (R,a)/2ir gives the degeneracy of the 
quasi-classical antidot level with given energy E = R 2 /2. 
The integral in Eq. ( |34| ) can be easily calculated, and it 
leads to the following number of degeneracy of quantum 
level corresponding to the classical two-circle trajectory: 

2E , 

N 2 (E, a) = — p 2 (a/V8E), (35a) 

P2 (x) = ((x)Q(l -x) + |e(.T - 1) - C(x/2), (35b) 

C(x) — arcsinx + x\/\ — x 2 . (35c) 



If the above number is less than unity, the level does not 
manifest itself. The energy values obtained by solving 
equation 

N 2 (E,a) = n, n= 1,2,3 (36) 

are indicated in Fig. || by dotted curves on its left side. 
We see that the larger is the energy the longer is the frac- 
tal plateau corresponding to the quasi-classical level, and 
the higher is its degeneracy. Moreover, there is a good 
coincidence of dotted curves with the plateaux ranges 
obtained by the numerical calculation. 

In a similar way the phase volume and the degeneracy 
of the quasi-classical levels corresponding to three circle 
classical trajectories can be estimated. In this case one 
can obtain 

n 3 (R, b) = e(2R - b)Q(a - 6)9(36 - a) 

x (a- 6)9(26 -a) + (36 - o)8(o - 26), (37a) 

N 3 (E,a) = ^p 3 (a/V8E), (37b) 

( 8((x/2) - 6((x/3) - 2C(as), < x < 1; 
p 3 (x) = I 8((x/2) - 6((x/3) -7T, 1< x < 2; 

[ 3tt- 6C(z/3), 2 < x < r. 

(37c) 

The energy values obtained by solving equation 

N 3 (E,a) = n, n = l,2 (38) 

are shown by dotted curves on the right side of Fig. || as 
well. We see that they also correlate well with the ranges 
of plateaux exceeding by 1/3 and 2/3 the corresponding 
Landau levels. 

A good agreement of all dotted curves with the degen- 
eracy of the fractional plateaux obtained in the numeri- 
cal solution of the problem convinces us of the adequacy 
of considered quasi-classical quantization scheme. The 
above picture is confirmed by the electron density plot 
presented in Fig. |j| The density is calculated for the 
third antidot level expelled from the fifth Landau state 
in the case ci/Ib — 6.32 indicated by small solid circle in 
Fig. 0. We see a rather good correlation of the electron 
density with the classical two-circle trajectory shown by 
a dotted curve in Fig. 0. 
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VII. CONCLUSIONS 



The energy spectrum of the electron moving in the per- 
pendicular magnetic field in a vicinity of impenetrable 
stripe and the corresponding densities are calculated by 
making use the integral equation technique for the anti- 
dot perimeter function (the perpendicular wave function 
derivative at the antidot border). It is shown that the 
perimeter function singularities caused by sharp edges 
of the antidot can be overcome by proper discretization 
technique which takes explicitly into account the loga- 
rithmic singularity of the kernel and root-type singulari- 
ties of the perimeter function. 

The antidot in the magnetic field expels the antidot 
energy levels from every degenerate Landau state. In the 
case of the round antidot due to the circular symmetry 
expelled antidot levels go up when the radius of the an- 
tidot increases (or the magnetic field strength increases) 
and freely intersect each other and the Landau levels. In 
the case of the antidot-stripe due to the lack of the sym- 
metry the variables can not be separated, and nearly all 
crossings are replaced by the anti-crossings. The expelled 
antidot levels can not cross the Landau levels, and conse- 
quently, they saturate below the next Landau level when 
the stripe length increases. 

As the antidot-stripe still has the inversion symmetry 
only the pairs of symmetric and anti-symmetric levels 
cross each other demonstrating the characteristic oscil- 
lations of the spectrum branches expelled from the ex- 
cited Landau levels. These oscillations and the above 
mentioned saturation can be explained in the asymptotic 
long stripe case by the interference of the magnetic edge 
modes. 

When the electron energy increases (for the antidot 
levels expelled from the higher Landau levels) the above 
mentioned oscillations of the spectrum branches is trans- 
formed into plateaux at the fractal cyclotron energy val- 
ues. It is shown that these fractal plateaux can be 
explained using simple quasi-classical quantization rule, 
and they are related to the classical trajectories com- 
posed of several Larmor circles. 

The above mentioned fractal plateaux of energy 
branches have to be seen in magnetoresitance of arrays 
of stripe type antidots, and the magnetization which is 
just proportional to the electron energy derivative over 
the magnetic field strength (in considered case through 
the dimensionless stripe length a). 



APPENDIX A: DISCRETIZATION OF 
SINGULAR INTEGRAL EQUATION 

In this section some details of the numerical solution 
of integral equation (0) are given. For the sake of conve- 
nience we scale the variables x ^ ax, and rewrite sepa- 
rately the equation for the symmetric and anti-symmetric 
perimeter function 



F ± (-x) = \{F(x)±F(-x)}. 
Now the equation reads 

f l/2 

/ dx'K ± (x,x')F ± (x')=0 
Jo 

with the symmetric (or anti-symmetric) kernel 
K^(x, x') — K(ax, ax') ± K(ax, —ax'). 



(Al) 



(A2) 



(A3) 



Next, we write down explicitly the singularity of the 
perimeter function 



F±(x) 



(A4) 



and the logarithmic singularity of the kernel 

K^fax') = 2\n\x-x'\+K ± {x,x') (A5) 

where the second part of the kernel K is regular func- 
tion at x = x' . It enables us to rewrite the integral 
equation as follows: 



1/2 rfx^lnlx-x'l/*^') 
^(1/2-30(1/2- a:') 
1/2 dx , A- ± (x,x , )/ ± (x ; ) 
o v /(l/2-x)(l/2- a; ') 



= 0. 



(A6) 



N ote both terms of the equation are divided by factor 
\/l/2 — x in order to have the symmetric kernel, as the 
numerical calculation of corresponding symmetric matrix 
eigenvalues can be performed with greater accuracy than 
those for the non symmetric one. 



The discretization of Eq. (A6) has been performed as 
follows. The perimeter function has been replaced by 
vector f* = {/o, /i, ■ • • } with components /„ = f ± (x n ), 
and x n = h(n + 1/2), h — 1/2N. Integral equation ( |Aq ) 
itself has been rewritten in the form of the following ma- 
trix equation: 
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{X+B ± }f ± = 0, 



(A7) 
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where the corresponding matrix elements of the kernel 
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are defined as 



An+l)h 

Anm / dx 

J 7ih 

X r 7 a^, (A8a) 

B± m = iPfa^BnBra, (A8b) 
r(n+X)h dx 



B n — 



n h 



(A8c) 



Both integrals can be calculated straightforwardly , and 
the analytical expressions for the discretization weights 



obtained. 

Using the discretized perimeter function f ± the elec- 
tron wave function (and the corresponding density) has 
been obtained via discretized version of Eq. (||). Calcu- 
lating the wave function outside the stripe (on the stripe 
it is equal to zero) only the perimeter function singular- 
ity has to be taken into account. Thus, Eq. (0) can be 
replaced by 



N-l 



^(r) = - B n {G(ar\ax',0)±G(ar\-ax',0)}f n . 

(A9) 



n=0 
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